Transfer Functions for Protein Signal Transduction: Application to a Model of Striatal Neural Plasticity

We present a novel formulation for biochemical reaction networks in the context of signal transduction. The model consists of input-output transfer functions, which are derived from differential equations, using stable equilibria. We select a set of 'source' species, which receive input signals. Signals are transmitted to all other species in the system (the 'target' species) with a specific delay and transmission strength. The delay is computed as the maximal reaction time until a stable equilibrium for the target species is reached, in the context of all other reactions in the system. The transmission strength is the concentration change of the target species. The computed input-output transfer functions can be stored in a matrix, fitted with parameters, and recalled to build discrete dynamical models. By separating reaction time and concentration we can greatly simplify the model, circumventing typical problems of complex dynamical systems. The transfer function transformation can be applied to mass-action kinetic models of signal transduction. The paper shows that this approach yields significant insight, while remaining an executable dynamical model for signal transduction. In particular we can deconstruct the complex system into local transfer functions between individual species. As an example, we examine modularity and signal integration using a published model of striatal neural plasticity. The modules that emerge correspond to a known biological distinction between calcium-dependent and cAMP-dependent pathways. We also found that overall interconnectedness depends on the magnitude of input, with high connectivity at low input and less connectivity at moderate to high input. This general result, which directly follows from the properties of individual transfer functions, contradicts notions of ubiquitous complexity by showing input-dependent signal transmission inactivation.Comment: 13 pages, 5 tables, 15 figure

Reconstructing the three-dimensional GABAergic microcircuit of the striatum

A system's wiring constrains its dynamics, yet modelling of neural structures often overlooks the specific networks formed by their neurons. We developed an approach for constructing anatomically realistic networks and reconstructed the GABAergic microcircuit formed by the medium spiny neurons (MSNs) and fast-spiking interneurons (FSIs) of the adult rat striatum. We grew dendrite and axon models for these neurons and extracted probabilities for the presence of these neurites as a function of distance from the soma. From these, we found the probabilities of intersection between the neurites of two neurons given their inter-somatic distance, and used these to construct three-dimensional striatal networks. The MSN dendrite models predicted that half of all dendritic spines are within 100 mu m of the soma. The constructed networks predict distributions of gap junctions between FSI dendrites, synaptic contacts between MSNs, and synaptic inputs from FSIs to MSNs that are consistent with current estimates. The models predict that to achieve this, FSIs should be at most 1% of the striatal population. They also show that the striatum is sparsely connected: FSI-MSN and MSN-MSN contacts respectively form 7% and 1.7% of all possible connections. The models predict two striking network properties: the dominant GABAergic input to a MSN arises from neurons with somas at the edge of its dendritic field; and FSIs are interconnected on two different spatial scales: locally by gap junctions and distally by synapses. We show that both properties influence striatal dynamics: the most potent inhibition of a MSN arises from a region of striatum at the edge of its dendritic field; and the combination of local gap junction and distal synaptic networks between FSIs sets a robust input-output regime for the MSN population. Our models thus intimately link striatal micro-anatomy to its dynamics, providing a biologically grounded platform for further study

The Nanostructure of Myoendothelial Junctions Contributes to Signal Rectification between Endothelial and Vascular Smooth Muscle Cells

Micro-anatomical structures in tissues have potential physiological effects. In arteries and arterioles smooth muscle cells and endothelial cells are separated by the internal elastic lamina, but the two cell layers often make contact through micro protrusions called myoendothelial junctions. Cross talk between the two cell layers is important in regulating blood pressure and flow. We have used a spatiotemporal mathematical model to investigate how the myoendothelial junctions affect the information flow between the two cell layers. The geometry of the model mimics the structure of the two cell types and the myoendothelial junction. The model is implemented as a 2D axi-symmetrical model and solved using the finite element method. We have simulated diffusion of Ca2+ and IP3 between the two cell types and we show that the micro-anatomical structure of the myoendothelial junction in itself may rectify a signal between the two cell layers. The rectification is caused by the asymmetrical structure of the myoendothelial junction. Because the head of the myoendothelial junction is separated from the cell it is attached to by a narrow neck region, a signal generated in the neighboring cell can easily drive a concentration change in the head of the myoendothelial protrusion. Subsequently the signal can be amplified in the head, and activate the entire cell. In contrast, a signal in the cell from which the myoendothelial junction originates will be attenuated and delayed in the neck region as it travels into the head of the myoendothelial junction and the neighboring cell

The Effects of NMDA Subunit Composition on Calcium Influx and Spike Timing-Dependent Plasticity in Striatal Medium Spiny Neurons

Calcium through NMDA receptors (NMDARs) is necessary for the long-term potentiation (LTP) of synaptic strength; however, NMDARs differ in several properties that can influence the amount of calcium influx into the spine. These properties, such as sensitivity to magnesium block and conductance decay kinetics, change the receptor's response to spike timing dependent plasticity (STDP) protocols, and thereby shape synaptic integration and information processing. This study investigates the role of GluN2 subunit differences on spine calcium concentration during several STDP protocols in a model of a striatal medium spiny projection neuron (MSPN). The multi-compartment, multi-channel model exhibits firing frequency, spike width, and latency to first spike similar to current clamp data from mouse dorsal striatum MSPN. We find that NMDAR-mediated calcium is dependent on GluN2 subunit type, action potential timing, duration of somatic depolarization, and number of action potentials. Furthermore, the model demonstrates that in MSPNs, GluN2A and GluN2B control which STDP intervals allow for substantial calcium elevation in spines. The model predicts that blocking GluN2B subunits would modulate the range of intervals that cause long term potentiation. We confirmed this prediction experimentally, demonstrating that blocking GluN2B in the striatum, narrows the range of STDP intervals that cause long term potentiation. This ability of the GluN2 subunit to modulate the shape of the STDP curve could underlie the role that GluN2 subunits play in learning and development

Soft-bound synaptic plasticity increases storage capacity

Accurate models of synaptic plasticity are essential to understand the adaptive properties of the nervous system and for realistic models of learning and memory. Experiments have shown that synaptic plasticity depends not only on pre- and post-synaptic activity patterns, but also on the strength of the connection itself. Namely, weaker synapses are more easily strengthened than already strong ones. This so called soft-bound plasticity automatically constrains the synaptic strengths. It is known that this has important consequences for the dynamics of plasticity and the synaptic weight distribution, but its impact on information storage is unknown. In this modeling study we introduce an information theoretic framework to analyse memory storage in an online learning setting. We show that soft-bound plasticity increases a variety of performance criteria by about 18% over hard-bound plasticity, and likely maximizes the storage capacity of synapses

Colocalization of Protein Kinase A with Adenylyl Cyclase Enhances Protein Kinase A Activity during Induction of Long-Lasting Long-Term-Potentiation

The ability of neurons to differentially respond to specific temporal and spatial input patterns underlies information storage in neural circuits. One means of achieving spatial specificity is to restrict signaling molecules to particular subcellular compartments using anchoring molecules such as A-Kinase Anchoring Proteins (AKAPs). Disruption of protein kinase A (PKA) anchoring to AKAPs impairs a PKA-dependent form of long term potentiation (LTP) in the hippocampus. To investigate the role of localized PKA signaling in LTP, we developed a stochastic reaction-diffusion model of the signaling pathways leading to PKA activation in CA1 pyramidal neurons. Simulations investigated whether the role of anchoring is to locate kinases near molecules that activate them, or near their target molecules. The results show that anchoring PKA with adenylyl cyclase (which produces cAMP that activates PKA) produces significantly greater PKA activity, and phosphorylation of both inhibitor-1 and AMPA receptor GluR1 subunit on S845, than when PKA is anchored apart from adenylyl cyclase. The spatial microdomain of cAMP was smaller than that of PKA suggesting that anchoring PKA near its source of cAMP is critical because inactivation by phosphodiesterase limits diffusion of cAMP. The prediction that the role of anchoring is to colocalize PKA near adenylyl cyclase was confirmed by experimentally rescuing the deficit in LTP produced by disruption of PKA anchoring using phosphodiesterase inhibitors. Additional experiments confirm the model prediction that disruption of anchoring impairs S845 phosphorylation produced by forskolin-induced synaptic potentiation. Collectively, these results show that locating PKA near adenylyl cyclase is a critical function of anchoring

Early phase of plasticity-related gene regulation and SRF dependent transcription in the hippocampus

Hippocampal organotypic cultures are a highly reliable in vitro model for studying neuroplasticity: in this paper, we analyze the early phase of the transcriptional response induced by a 20 \ub5M gabazine treatment (GabT), a GABA-Ar antagonist, by using Affymetrix oligonucleotide microarray, RT-PCR based time-course and chromatin-immuno-precipitation. The transcriptome profiling revealed that the pool of genes up-regulated by GabT, besides being strongly related to the regulation of growth and synaptic transmission, is also endowed with neuro-protective and pro-survival properties. By using RT-PCR, we quantified a time-course of the transient expression for 33 of the highest up-regulated genes, with an average sampling rate of 10 minutes and covering the time interval [10 3690] minutes. The cluster analysis of the time-course disclosed the existence of three different dynamical patterns, one of which proved, in a statistical analysis based on results from previous works, to be significantly related with SRF-dependent regulation (p-value<0.05). The chromatin immunoprecipitation (chip) assay confirmed the rich presence of working CArG boxes in the genes belonging to the latter dynamical pattern and therefore validated the statistical analysis. Furthermore, an in silico analysis of the promoters revealed the presence of additional conserved CArG boxes upstream of the genes Nr4a1 and Rgs2. The chip assay confirmed a significant SRF signal in the Nr4a1 CArG box but not in the Rgs2 CArG box

25th Annual Computational Neuroscience Meeting: CNS-2016

Abstracts of the 25th Annual Computational Neuroscience Meeting: CNS-2016 Seogwipo City, Jeju-do, South Korea. 2–7 July 201